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ABSTRACT 

We present the analysis of Spitzer-IRS spectra of four early-type galaxies, NGC 1297, NGC 5044, NGC 6868, and 
NGC 7079, all classified as LINERs in the optical bands. Their IRS spectra present the full series of H2 rotational 
emission lines in the range 5-38 /jm, atomic lines, and prominent PAH features. We investigate the nature and origin 
of the PAH emission, characterized by unusually low 6 — 9/11.3 /im inter-band ratios. After the subtraction of a 
passive early type galaxy template, we find that the 7 — 9 /im spectral region requires dust features not normally 
present in star forming galaxies. Each spectrum is then analyzed with the aim of identifying their components and 
origin. In contrast to normal star forming galaxies, where cationic PAH emission prevails, our 6-14 /im spectra seem 
to be dominated by large and neutral PAH emission, responsible for the low 6 — 9/11.3 /im ratios, plus two broad 
dust emission features peaking at 8.2 /im and 12 /im. Theses broad components, observed until now mainly in evolved 
carbon stars and usually attributed to pristine material, contribute approximately 30-50% of the total PAH flux in 
the 6-14 /im region. We propose that the PAH molecules in our ETGs arise from fresh carbonaceous material which 
is continuously released by a population of carbon stars, formed in a rejuvenation episode which occurred within the 
last few Gyr. The analysis of the MIR spectra allows us to infer that, in order to maintain the peculiar size and charge 
distributions biased to large and neutral PAHs, this material must be shocked, and excited by the weak UV interstellar 
radiation field of our ETGs. 

Subject headings: Galaxies: elliptical — Galaxies: evolution — Galaxies: stellar content — Galaxies: 
individual (NGC 1297, NGC 5044, NGC 6868, NGC 7079) 



1. INTRODUCTION 

The understanding of the evolution of early-type galax- 
ies (i.e. E and SO galaxies) has been greatly enhanced 
by the study of their interstellar medium (ISM) (see e.g. 
Rcnzini 2007, for a review). This component, and its 
relevance in early-type galaxies (ETGs hereafter) was 
widely neglected in early studies, since these galaxies 
were, for a long time, considered to be essentially de- 
voided of interstellar gas and dust. In the last two 
decades, however, multi- wavelength observations have 
changed this picture and have detected the presence of a 
multi-phase ISM (e.g. IBregman et al] 1992: B ettoni et al.l 
l200lUBoselli et al.f 2005: M organti et al.ll2006[ ). The bulk 
of the gas in ellipticals is h eated to the virial tem perature, 
emitting in X-rays (e.g. iBregman et al.l Il992[ ). Some 
ETGs also contain large dust masses, which is surprising 
given the efficiency of sputtering by their hot plasma (e.g . 
Knapp et al.lfl9 89: Ferr ari et al.lll999b iTemi et al.l 120041 : 
Clemens et all 12010ft . Even more surprising was the 
detection of polycycl ic aromatic hydrocarbons (PAHs) 
in ETG spectra (e.g. IXilouris et all 120041: iKaneda et all 
120051 120081 : IBressan et al.l l2006at iPanuzzo et al.l 120071 : 
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IBregman et al.l 120081) . because, dust molecules as small 
as PAHs, are expected to be more eas ily destroyed by 
the in teraction with the hot plasma (e.g. lDwek fc Arendtl 
fl992l) . 

PAH emission is a typical feature of the mid- 
infra red (MIR hereafter) s pectra of star-forming ga laxies 
(e.g. iRoussel et all l200lt iForster Schreiber et all 12004 
iPeeters et al.l l2004t I Vega et ahl 120051 ) , which are domi- 
nated by the strong 6.2, 7.7, and 8.6 /im emission bands 
usually attributed to the C-C stretching and C-H in- 
plane bending modes of ionized PAHs (see Tielens 2008 
and references therein). However, the PAH spectra of 
ETGs are often different from those of star forming 
galaxies. The most striking characteristic is the un- 
usual PAH inter-band emission ratio between the 7.7 //m 
complex and the 11. 3 /im complex (jKaneda et al.l 120051 : 
IBressan et al.l l2006a| ) . So, while in star forming galax- 
ies the 6.2, 7.7, and 8.6 /im features dominate the PAH 
spectrum, they are very weak in these ETGs, in spite of 
a strong 11.3 /im band. 

There already exist studies in the literature aimed 
at understanding the nature and origin of these PAH 
emission features, as well as the physical conditions of 
the ISM, that allow the surviva l of th is material in 
the observed ratios. ISmith et al.l (|2007) analyzing the 
SINGS galaxies sample, found that some LINERs also 
have unusual PAH ratios. They suggest that these ra- 
tios could be due to the relatively hard radiation field 
from a low-luminosity AGN, which modifies the emit- 
ting gr ain distribution at th e expense of the smaller 
PAHs. Breg man et al.l (|2008| ) claim that, once the con- 
tinuum of the underlying stellar population, which has 
a broad dip just below the 7.7/zm complex, is properly 
subtracted, the low 7.7 to 11.3 /im ratio turns into a 



2 



Vega et al. 



value typical of the I SM of late-type galaxies. In contrast, 
iKaneda et al.l (|2008| ) have shown cases where these ratios 
remain anomalous even after the proper removal of the 
underlying stellar population. These authors conclude 
that the anomalous ratios are likely due to the presence 
of a larger fraction of neutral PAHs with respect to ion- 
ized PAHs, contrary to what happens in the ISM of late 
type galaxies. They ascribe this anomalous mixture to 
recent accretion events. 

In this paper we deepen the investigation of this 
problem by presenting a new approach to the study 
of these peculiar ETGs. We analyze our, high S/N 
Spitzer-IRS observations of four ETGs (see Table IA7I 
for the galaxy classification) in low density environments 
(LDEs), namely NGC 1297, NGC 5044, NGC 6868, and 
NGC 7079, which display "unusual" PAH emission in 
their MIR spectra. For our study we make use of theo- 
retical spectra and empirical dust feature templates re- 
cently reported in the literature (see e.g. Rapacioli et al. 
2005, Berne et al. 2007, Joblin et al. 2008; Bauschlicher 
et al. 2008, Joblin et al. 2009, Berne et al. 2009). 
This method has the advantage of isolating the compo- 
nents responsible for the emission, thus providing more 
insight i nto the origin and evol ution of this ISM com- 
ponent (|Allamandola et ail |1999() . We advance a new 
scenario, where this PAH emission is naturally explained 
in terms of the presence of an underlying intermediate 
age population of carbon stars, which are continuously 
feeding the ISM with pristine carbonaceous material. If 
our interpretation is correct, the anomalous PAH spec- 
tra are strong evidence of a rejuvenation episode that 
affected the evolution of the ETGs in the recent past. 

The paper is organized as follows. In § 2 we provide the 
details of the observations and data reduction methods. 
In § 3 we measure and provide detailed information for 
each single PAH emission feature observed in the spectra 
and for PAH inter-band ratios. This information is used 
to compare our spectra with existing classifications. In 
§ 4 we explore the nature of the PAH carriers by using 
component templates available both from observations 
and theory. This method allow us to identify families 
of emitters typical of different environments. In § 5 we 
analyze the MIR emission lines sensitive to the physical 
conditions of the medium giving place to the PAH emis- 
sion. In § 6 we discuss our results and compare them with 
other studies. The relevant properties of our ETGs, as 
known so far, are provide in the Appendix. 

2. MID INFRARED OBSERVATIONS AND DATA 
REDUCTION 

Spitzer IRS spectral observations of all galaxies were 
obtained during the third Spitzer General Observer Cy- 
cle 3 on 2007 June 1 as part of program [ID 30256; P.I.: 
R.Rampazzo]. 

The observations were performed in Standard Staring 
mode with low resolution (R ~ 64-128) modules SL1 
(7.4-14.5/xm), SL2 (5-8. 7 urn), LL2 (14.l-21.3pm), and 
LL1 (19.5-38^m). Exposure times for each galaxy are 
given in Table []] 

We have devised an ad hoc procedure to flux calibrate 
the spectra that exploits the large degree of symme- 
try that characterizes the light distribution in early-type 
galaxies. The reduct ion procedure is fully described in 
iBressan et al.l (|2006al) ; here we only summarize the main 



TABLE 1 

IRS OBSERVATIONS 



Galaxy 
IRS module 


ramp duration number of cycles 
[seconds] 


AOR Key 
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Q 1 "1 


60 


19 




SL2 


60 


19 
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120 


14 




LL1 


120 


8 
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60 


19 
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60 


19 
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120 


11 
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120 
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CT 1 


60 


6 
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60 


6 
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120 


13 




LL1 


120 


8 




NGC 7079 
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SL1 


60 


19 




SL2 


60 


19 




LL2 


120 


11 




LL1 


120 


8 





steps. 

For SL observations, the sky background was removed 
by subtracting observations taken in different orders, 
but at the same nod position. LL segments were sky- 
subtracted by differencing the two nod positions. 

We obtained new e _ s _1 to Jy flux conversions by 
applying a correction for aperture losses (ALCF) and 
a correction for slit losses (SLCF) to the flux conver- 
sion tables prov i ded b y the Spitzer Science Center (e.g. 
iKennicutt et al.l I2003D . By applying the ALCF and 
SLCF corrections we obtained the flux received by the 
slit. 

For each galaxy we then simulated the corresponding 
observed profile along the slits by convolving a wave- 
length dependent bidimensional intrinsic surface bright- 
ness profile with the instrumental point spread function 
(PSF). The a dopted profile was a two dimensional modi- 
fied King law (jElson et al.lll987l ). By fitting the observed 
profiles with the simulated ones, we could reconstruct the 
intrinsic profiles and the corresponding intrinsic spectral 
energy distribution (SED). This procedure allowed us to 
determine whether a particular feature was spatially ex- 
tended or not. 

Finally the spectrum was extracted in a fixed width 
(18" for SL and 10.2" for LL) around the maximum in- 
tensity. After normalization to the same area (3.6" x 18") 
there remain differences between the SL and LL spec- 
tral segments of between 3% (NGC 5044) and 14% 
(NGC 6868). In order to account for these differences, 
the LL spectrum has been rescaled to match the SL one 
in the common wavelengths. 
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Fig. 1— Spitzer-YRS spectra of NGC 1297, NGC 5044, 
NGC 6868, and NGC 7079. Dotted lines indicate the positions 
of the most relevant PAH complexes at 6.2 fim, 7.7 /im, 8.6 (im, 
11.3 /im, 12.7 /im, and 17 /im. Arrows and thin solid lines indicate 
the positions of the H2 rotational lines. For comparison, we also 
show the spectra of an ETG with a recent burst of star formation 
(NGC 4435, Panuzzo et al. 2007) and of two passively evolving 
ETGs: the cluster ETG, NGC 4621 (Bressan et al. 2006a), and 
the field ETG, NGC 3818 JPanuzzo et al. l201Ch . both scaled to the 
flux of NGC 1297 at 5.6/iin. 



TABLE 2 

2MASS MAGNITUDES WITHIN THE CENTRAL 5 ARCSEC RADIUS 



NAME 


J 


H 


K 


NGC 1297 
NGC 5044 
NGC 6868 
NGC 7079 


11.98±0.01 
11.24±0.01 
10.41±0.01 
ll.10i0.01 


11.24±0.01 
10.51±0.01 
9.67±0.01 
10.41±0.01 


11.Q1±0.Q1 
10.22±0.01 
9.40±0.01 
10.17±0.Q1 



The uncertainty on the flux was evaluated by consid- 
ering two sources of noise: the instrumental plus back- 
ground noise and the Poissonian noise of the source. The 
former was evaluated by measuring the variance of pixel 
values in background-subtracted, co-added images far 
from the source. The Poissonian noise of the sources 
was estimated as the square root of the ratio between 
the variance of the number of e~ extracted per pixel in 
each exposure, and the number of the exposures. The 
total noise was obtained by summing the two sources in 
quadrature and by multiplying by the square root of the 
extraction width in pixels. We notice that the overall 
absolute photometric uncertainty of IRS is 10 per cent, 
while the slope deviation within a single segment (affect- 
ing all spectra in the same way) is less than 3 per cent 
(see the Spitzer Observers Manual). 

The flux calibrated de-redshifted IRS spectra of 
NGC 1297, NGC 5044, NGC 6868, and NCG 7079 are 
shown in Figure [JJ In the same figure we also plot, 
for comparison, the spectra of an ETG which has ex- 
perienced a recent ep isode of star formation, NGC 4435 
(jPanuzzo et al.ll2007| ) , and of two typical p assively evolv- 
ing ET Gs, the cluster ETG NGC 4621 (se elBressan et al. 
2006al) and the field ETG NGC 3818 (jPanuzzo et al. 
2010h . In spite of the similarity to the passive galax- 
ies, especially in the lower wavelength spectral range, 
the IRS spectra of our objects show several signatures of 
activity, such as: (a) strong PAH emission at Aq ~ 11.3 



/im, similar to that found bv lKaneda et al.l ()2005l ) in four 
ETGs and the well developed complex of PAH emission 
near 17 /im; (b) other weak dust features at 6.2 /im, 7.7 
jum, 8.6 /im, a nd 12.7 /im, usually att r ibuted to PAH 
molec ules (e.g . lAllamandola et aD 119891 : iPuget fc Legerl 
119891 : iTielensI 120081 ) ; (c) strikingly strong emission of 
pure H2 rotational lines; (d) forbidden nebular emission 
lines of Ar, Fe, Ne, O, S, and Si; and (e) warm dust 
continuum emission. 

3. ANALYSIS OF THE PAH EMISSION FEATURES 

3.1. Removal of the underlying stellar and dust continua 

In order to study the PAH features it is necessary to 
properly remove from each spectrum the underlying stel- 
lar and the dust continua. 

We adopt a semi-empirical, high S/N, template derived 
from the analysis of passive ETGs to describe the stel- 
lar continuum representative of the old stellar popula- 
tion. That template was built averaging the NIR (J-H-K 
2MASS) data and the 5-40 /im IRS-Spitzer spectra of 
three passively evolving ETGs in low density environ- 
ments, namely NGC 1389, NGC 1426, and NGC 3818 
(jPanuzzo et al.l [20101 ). The stellar continuum template 
shows the dip at 8 /i m likely due to photo spheric SiO ab- 
sorption bands (e.g. lVerhoelst et al.ll2009l ) and the bump 
of silicate emission at «10 /im likely due to dusty circum- 
stellar O-rich AGB envelop es. Both are feature s typica l 
of passively evolving ETGs (jBressan et alj|1998ll2006al) . 
In Figure [TJ as an example, we illustrate the range of 
variation of these stellar features in the stellar contin- 
uum of NGC 4621 and NGC 3818 located in Virgo and 
in the field, respectively. 

In order to properly remove the underlying continuum 
from the IRS spectra, we use the 2MASS magnitudes 
within an aperture of 5" radius, which is close to that 
of our IRS spectra dPanuzzo et al.ll2007l ). The NIR mag- 
nitudes and errors are listed in Table [3J Assuming that 
the NIR fluxes are completely due to the stellar compo- 
nent, we normalize our stellar continuum template to the 
observed flux in the H-band, and calculate the contribu- 
tion of the stellar continuum to the MIR spectra. We 
find the old stellar population contributes about 97.5%, 
98.5%, 99.6%, and 98.5% of the observed fluxes at 5.6 /mi 
for NGC 1297, NGC 5044, NGC 6868, and NGC 7079, 
respectively. We then may consider that the old stel- 
lar population dominates the MIR emission at shorter 
wavelengths. This was already expected from the com- 
parison with passive ETGs shown in Fig. [TJ However, 
this is not always the case. Indeed Panuzzo et al. (2007) 
found that the contribution of the stellar component to 
the MIR emission of NGC 4435 was less than 30% at 
5.6 /im. 

After a proper subtraction of the stellar continuum, 
the galaxies show a residual continuum component due 
to dust emission. It is removed by fitting up to 5 modi- 
fied black-bodies, with temperatures ranging from 20 to 
350 K. Notice that these dust components are chosen 
just to smoothly reproduce a realistic underlying dust 
continuum without any attempt to provide a physical 
interpretation. 

The whole process of continuum subtraction and the 
final results are shown in Figure [2] The final spectra are 
a combination of emission lines and dust features. 

All the spectra, normalized to the peak of the 11.3 /im 
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TABLE 3 

Parameters and integrated fluxes of the fitted PAH features. 
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Fig. 2. — Detailed fit to the MIR continuum emission of our sample of ETGs. The orange line represents the template of the passive 
ETG normalized to the H-band. Blue lines represent the different thermal dust components. The green line represents the sum of all dust 
components. The red line is the best fit model to the total continuum emission, calculated as the sum of the normalized passive ETG plus 
the thermal dust components. In the insets we plot the continuum subtracted spectra. Notice the broad feature around 8 /im that emerges 
when the continuum is subtracted. 
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Fig. 3. — Continuum-subtracted spectra of our ETGs normal ized to the peak int ensity of the 11.3 fJ.va feature. For comparison, we also 
show the IRS spectrum of NGC 1482, a SINGs HII galaxy from Smith ct al. (2007). Vertical lines indicate the positions of the main PAH 
complexes at 6.2, 7.7, 8.6, 11.3 12.8, and 17 fim. 



feature, are shown in Figure [3] For comparison, we also 
plot in the same figure the spectrum of NGC 1482, a 
HII gal axy with PAH emi ssion typical of a star-forming 
galaxy ([Smith et al.ll2007tl . 

3.2. Analysis of the dust feature emission 

Before discussing the individual PAH features, we 
point out the following points on the general shape of 
the spectra: 



• The positions of the peak of the 7.7 /mi band in 
NGC 6868 and in NGC 5044, are shifted to longer 
wavelengths. 

• Despite the similarity of the PAH emission bands 
in the 4 ETGs, the emission lines show consider- 
able variation from galaxy to galaxy. These emis- 
sion lines will be analyzed for a larger sample in a 
companion paper ([Panuzzo et al.ll2010t ). 



The continuum subtracted spectra of our ETG 
sample are very similar. At the same time, in the 5— 
9 /mi spectral range, they significantly differ from 
star forming galaxies, represented by NGC 1482. 

A broad 7.7 /im feature, hardly recognizable in the 
original spectra, appears after continuum subtrac- 
tion because of the presence of the dip at 8 /mi in 
the underlying stellar continuum. 

The PAH bands at 8.6, 7.7, and primarily at 6.2 /mi 
are much less intense than those observed in the 
integrated spectrum of star forming galaxies (see 
the NGC 1482 spectrum). On the contrary, the 
emission at 17 /tm looks very similar to that of star 
forming galaxies. 



3.2.1. The individual features 

We use a spe ctral decomposition method similar to 
that outlined by ISmith et all ([2007( 1 to recover the PAH 
features and the emission lines. In this method, the 
spectrum is modelled by a combination of Gaussian pro- 
files for the atomic and molecular emissio n lines, and 
Prud e profiles for the PAH features (e.g. iLi fc Draind 
12001( 1. The broad PAH features are represented by up to 
four sub-features grouped as blended complexes, some 
of them adopted to best reproduce the band shape (e.g. 
Smith et al. 2007). In order to explore if the peculiarities 
of our PAH spectra could be due not only to differences 
in the relative intensities of the sub-features, but also to 
differences in their characterizations, we leave as free pa- 
rameters the central wavelengths, An and the fractional 
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Fig. 4. — Detailed fit to the lines and the PAH emission for the continuum subtracted spectra of NGC 1297 (upper left), NGC 5044 
(upper right), NGC 6868 (lower left), and NGC 7079 (lower right) in the spectral region between 5 and 20 fim. The open squares and the 
blue solid line indicate the observed MIR continuum subtracted spectra. The solid red line indicates the total fit to the data, calculated as 
the sum of the PAH features (solid thin green lines), and the emission lines (solid thin orange lines). Arrows indicate the positions of the 
fitted PAH features, nebular emission lines, and H2 rotational lines. 
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TABLE 4 

PAH INTERBAND RATIOS AND PAH TO CONTINUUM RATIOS 



Lpah(6.2)/«/Xv(24) 
Lpah(7.7)/i/£v(24) 
L P ah(11.3)/^L„(24) 
L pah (17)/i/L„(24) 



N1297 



N5044 



N6868 



N7079 



0.20 ±0.11 
0.27 ±0.10 
0.54 ±0.31 
1.36 ±0.09 
3.40 ± 0.36 
2.72 ±0.20 
0.31 ±0.04 
2.19 ±0.12 
0.50 ±0.01 



0.15 ± 0.06 
0.47 ±0.11 
0.51 ± 0.21 
1.80 ±0.13 
3.77 ±0.28 
3.34 ±0.26 
0.47 ±0.03 
2.18 ± 0.17 
0.44 ± 0.02 



0.18 ±0.04 
0.40 ± 0.09 
0.77 ±0.15 
2.18 ±0.09 
4.14 ±0.23 
4.18 ±0.18 
0.53 ±0.02 
1.90 ±0.13 
0.52 ±0.01 



0.23 ±0.03 
0.30 ±0.04 
0.87 ±0.11 
1.26 ±0.10 
4.79 ±0.55 
3.78 ±0.33 
0.27 ±0.02 
2.03 ±0.14 
0.35 ±0.02 



0.058±0.034 0.039±0.009 0.044±0.010 0.031±0.005 

0.220±0.029 0.148±0.010 0.240±0.025 0.132±0.017 

0.215±0.022 0.082±0.008 0.110±0.011 0.101±0.010 

0.107±0.011 0.038±0.004 0.057±0.006 0.035±0.004 



SINGs median valuc^ 



0.28 
1.10 
1.90 
3.60 
5.70 
6.90 
0.68 
1.80 
0.53 



0.10 
0.40 
0.10 
0.06 



0.21 - 0.34] 
0.77 - 1.43] 
1.33 - 2.47] 
2.52 - 4.68] 
3.99 - 7.41] 
4.83 - 8.97] 
0.48 - 0.88] 
1.26 - 2.39] 
0.37 - 0.69] 



0.05 - 0.20] 
0.25 - 1.00] 
0.04 - 0.20] 
0.03 - 0.10] 



Notes: a blended PAH complex, see Table [3] ^ 



median value and 1 <r variation in the corresponding ratio derived by Smith ct al. (2007) 
for the SINGs galaxy sample. 
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full width at half maxima (FWHM), 7, of the Drude pro- 
files. The model is then fitted to the observed spectra by 
minimizing the global y 2 using the Levenberg-Marquardt 
algorithm, and usi ng the values Ao and 7 reported by 
iDraine fc Lil (|2007t ) (hereafter DL07) as the initial guess 
for the parameters of the Drude profile. 

The best fits to the spectra are shown as thick red 
lines in Figure HI where the fitted PAH subfeatures and 
emission lines are shown as thin green and thin orange 
lines, respectively. The results of the fits are summarized 
in Table [3] Columns 1 and 2 list the central wavelengths, 
A and the fractional FWHM, 7, reported by DL07. 

In the remaining columns we report, for each galaxy, 
the fitted Ao, 7. and the corresponding integrated fluxes. 
Different sub-features of the 7.7 /im, 11.3 /mi, 12.7 /im, 
and 17 /jm PAH complexes, are indicated with the let- 
ters a, b, c, and d, respectively. The errors in the derived 
fluxes and Drude profile parameters have been estimated 
repeating the above procedure with simulated spectra 
obtained by randomly varying the fluxes within the ob- 
servational error. The process was repeated 100 times 
and the errors were calculated as the standard deviation 
of the distribution of the fitted values in the simulations. 

From Table [3] we notice that even though most of the 
Ao and 7 of the fitted features are in good agreement 
with those reported by DL07, the fits for NGC 6868, 
NGC 5044, and to a lesser extent NGC 1297, required 
some Drude profiles with significantly different lambda 
peaks and/or fractional FWHM. For some features the 
differences may be due to low signal to noise (i.e. 5.27 
/im) but this is not the case for the 7-9 /im complex, 
which is well detected in all cases. 

The most remarkable differences with respect to DL07 
are: 

• The weakness or even lack of the 7.6 /im feature 
in NGC 5044 and NGC 6868. The 7.6 /im feature 
is the main component of the 7.7 /im complex in 
normal and active galaxies (e.g. iPeeters et al.ll2004j : 
iSmith et al.H2007t rTielensll2008l) ~ 

• The brightest component is the 7.85 /im sub- 
feature. It is also about twice as wide (7 ~ 
0.077) as that in DL07. Moreover, in the case of 
NGC 5044, it peaks at A = 7.90/tm. 

• The fits to the broad 6-9 /im feature in the 
NGC 5044 and NGC 6868 spectra require also of 2 
times wider (7 ~ 0.75 — 0.8) Drude profiles for the 
8.6 /im feature. 

• The 17 /im PAH complex is also wider than in 
normal galaxies and requires an "extra" feature at 
about 19.1 /im. 

The fitting requirements of wider Drude profiles for 
the well established features at 7.7 /im and 8.6 /tm are 
actually due to an emission excess in the 6-9 /im spec- 
tral range with respect to what is recovered by using 
only the standard PAH Drude profiles (i.e. DL07's fea- 
tures). Thus, the broad 6-9 /im feature in NGC 5044 
and NGC 6868 could also be well reproduced by adding 
some extra Drude profile sub-features to the standard 
PAH feature set. 

In order to test if this excess emission could be due 
to the use of a different decomposition method to that 



used by Smith et al. (2007), we ap plied our procedure 
to NGC 4435 (iPanuzzo et all 120071) and to NGC 1482 
(|Smith et a l. 2007) from the Smith et al. sample. In nei- 
ther case did we find significant differences with respect 
to the values reported by DL07. By using the NGC 1482 
spectrum we also test the dependence of our results on 
the different methods of subtraction of the underlying 
continuum. Following the recipe described in Section 3 
we find that the contribution from the old stellar popula- 
tion to the 5.6 /tm flux is < 28%, and that the differences 
between the continuum subtracted spectrum obtained 
with our method and with PAHFIT (Smith et al. 2007) 
are negligible. Thus, we are confident that our results are 
not an artifact of the adopted method for the subtraction 
of the underlying continuum. Notice however, that for 
a correct subtraction of the stellar component, the use 
of a proper stellar template is mandatory for ETGs at 
odds with HII galaxies where the contribution of the old 
stellar population to the MIR spectral range is small. 

3.2.2. PAH inter-band ratios 

The strengths of the individual PAH bands are affected 
by a variety of mechanisms, that include the distribution 
of PAH sizes and charges, PAH de-hydrogenization, and 
the strength a nd hardness of the radiation field e xciting 
the PAHs (e.g. iTielens! 120081: iGalliano et al.ll2008f ). Pos- 
sible destruction mechanisms of PAHs are also varied, in- 
cluding pho to-dissociation, sputtering and dissociation in 
shocks (e.g. iLeger et al.lll989t iMicelotta et al.ll2010allbT) . 
Variations of the PAH inter-band ratios can be explained 
by the combined effects of those mechanisms, and hence 
to variatio ns in the physical conditions of the environ- 
ment (e.g. lAllamandola et al.lll989t ISchutte et al.l 119931: 
IDraine fc Lill2007tlGalliano et ; al.ll2008l : lTielensll2008h . 

Table U lists the ratios of the strongest PAH bands. 
In the last column of the table we report the median 
values and their 1 a va riation for the 59 o bjects of the 
SINGs galaxy sample (|Smith et al.l l2007| ). The most 
striking characteristic of the spectra is the unusually low 
7.7/11.3 /tm PAH inter-band ratio with respect to what 
is observed in star forming g alaxies. This unusu al ra- 
tio was already repo rted by iKaneda et al.l (|2005l ) and 
iBressan et al.l (2006a) i n a small set of ETGs and later 
confirmed by|Kaneda et al. (2008) for a larger sample of 
ETGs and by ISmith et al.l (|2007| ) for a subset of about 
10 objects in the SINGs galaxy sample. A more careful 
inspection of the spectra indicates that the 7.7/11.3 /tm 
PAH inter-band ratio is not alone in presenting these 
unusual values but there is a systematic deficit in the 
shorter wavelength bands (i.e. 6.2, 7.7, and 8.6 /im) 
with respect to those peaking at longer wavelengths (i.e. 
11.3, 12.7, and 17 /tm). This effect is more extreme at 
shorter wavelengths, with the 6.2 /tm band having the 
strongest deficit. The average 6.2/11.3 ratio for our four 
galaxies is 0.36 ±0.20 which is significantly different from 
1.10 ± 0.33 of the SINGs galaxies. The same for 6.2/17, 
our is 0.67 ±0.30 while SINGs is 1.9 ±0.57. Differences in 
the 7.7/11.3 and 7.7/17 /tm ratios are less extreme. Our 
average 7.7/11.3 and 7.7/17 /tm ratios are 1.65±0.54 and 
3.51 ± 0.80, while the mean values for the SINGs sam- 
ple are 3.60 ± 1.08 and 6.90 ± 2.07, respectively. In the 
case of the 8.6 /im band, only NGC 1297 and NGC 7079 
present 8.6/11.3 /tm and 8.6/11.3 /tm ratios which are 
significatively lower than the values of the SINGs sample. 
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Whereas, 17/11.3 /im and the 12.7/11.3 /im ratios are, 
in all the cases, within the values found for the SINGs 
galaxy sample. We check wether these differences could 
be due to some general PAH emission defect with respect 
to dust emission. In the last rows of Table|4]we report the 
ratio between PAH complexes and the MIPS24 / j m lumi- 
nosity. The last column also provides the median values 
of the SINGs sample. We see again that the deficit is 
selective, and it occurs only in the shorter wavelength 
bands, while at 11.3/im and 17/im, the ratios fall within 
those of the SINGs sample. 

Laboratory studies and theoretical calculations show 
that the 6 — 9/im bands are intrinsically weak in neu- 
tral PAHs and become stronger when the PAHs are ion- 
ized (Allamandola et al. 1999; Bauschlicher 2002; Kim 
& Saykally 2002). Therefore, the 6 — 9/im bands will be 
much more intense for an ionized PAH than for a neutral 
PAH, while the contrary will be true for the 3.3 and 11.3 
/xm bands. Consequently, the ratios between the 6 — 9/im 
bands and the 11.3 /jm band depend on the charge of the 
PAHs, which is directly related to the physical conditions 
(intensity of the ionizing radiation field, electron density, 
etc.) of the PAH emitting region. On the other hand, 
larger PAHs attain lower temperatures than their smaller 
counterparts upon the absorption of the same photon. 
That implies that smaller PAHs will contribute mainly 
at short wavelengths, while larger PAHs a ccount mainly 
for the emi s sion at long wavelengths (e.g. ISchutte et all 
19901 119931 iDrainefcjijjgOp: iB auschlic her et al I I2008L 
2009UBoersma et aJJ[2010h . iBoersma et all (|2010t l found 
that while the 6 — 9 /im and 11.3 /im bands are very sen- 
sitive to the ionization state of the PAH molecule, the 17 
/im complex is not, and that the emission of this com- 
plex is probably due to the C-C-C vibrational modes of 
a limited number of large PAHs. Consistently, the ratios 
of the 17/6 — 9 /im features would be a sensitive probe 
of the size of interstellar PAHs so that size distributions 
biased to large PAHs would result in higher values of the 
17/6.2 /im and 17/7.7 /im ratios compared to the values 
in normal galaxies. On the other hand, charge distribu- 
tions biased to neutral PAHs will result in lower 7.7/11.3 
/im and 6.2/11.3 /im ratios compared to those in normal 
galaxies. 

Therefore, the low values of the 6.2/17 /im and the 
7.7/17 /im ratios in our ETGs imply a size distribution 
deficient in small PAHs, and responsible for the "nor- 
mal" PAH emission at 17 /im. This suggestion is rein- 
forced by the small values of the 6.2/7.7 /im ratio which 
is particularly sensiti ve to the lower cutoff of the PAH 
size d istribution (e.g. iDraine "fcH 120071 iGalliano et all 
2008). However, theoretical computations of large ion- 
ized PAH spectra indicate an enhanced emission of the 
6 — 9/im bands with respect to the 11.3 /im band, even 
if the effects temperature cascade is tak en into account 
in th eir computations (see Figure 17 of IBoersma et all 
l2010h . Therefore, a mixture of large PAHs with a normal 
charge distribution could not reproduce the low values 
of the 7.7/11.3 /im ratios observed in our objects, which 
would require of a mixture plentiful of neutral PAHs. In 
summary, the low values of the 6 — 9/im bands versus the 
11.3 /im and the 17 /im bands observed in our objects 
are compatible with the emission from a mixture which 
is rich in large and neutral PAHs. 



4. CHARACTERIZATION OF THE PAH EMISSION 

Instead of using combinations of single Drude pro- 
files, we will now analyze the continuum subtracted spec- 
tra by means of emission templates of different dust 
components. We wi ll follow the metho d devised by 
IRapacioli et al.l (I2005D . iBerne et al.l (120071) [Joblin et all 
(|2008l ) . Uoblin et alJ(|2009f ) and lBerne et alJ(|2009t l which 
provides general information on the sizes and ionization 
states of PAHs, and the presence of other dust compo- 
nents. In this way, we will gain insight into the nature 
of these spectral features and, ultimately, about the stel- 
lar populations and physical processes involved. Before 
using the emission templates, we removed the emission 
lines from the continuum subtracted spectra by using the 
values of the lines derived in the fits performed in section 



4.1. The Templates 

There are several studies in the literature about emis- 
sion templates of PAH families and other dust features 
(e.g. ICesarskv et all 120001: IRapacioli et al l 1 20051. 120061: 
Berne et all l2007t Uoblin et all I2008L 12009b Berne et all 
2009h . For the following ana lysis we will u se the PAH- 
dust feature templates from Ijoblin et all (|2008l ) which 
are shown in Figure [5] The parameters of the tem- 
plates (i.e. Ao and FWHM of the Lorentzian and Gaus- 
sian profiles) were taken from their Table A.l. Three 
of them were obtained from the mathematical decom- 
position of three reflection nebulae observed with ISO 
and Spitzer b y using the Principal Component Analy - 
sis technique (jRapacioli et al.l 120051: IBerne et all [20071 
The first template, called PAH + , characterized by strong 
C-C modes (6-9 /im) relative to C-H modes (10 - 
14 /im) and peaking at 7.6 /i m, is attributed to PAH 
cations (Rap acioli et al.l 120051) . The second template 
called PAH , characterized by strong C-H modes relative 
to C-C modes and showing an evident peak at 11.3 /im, is 
attributed by the authors to a PA H mixture rich in neu- 
tral PAHs (|Rapacioli et all 1200511 The third one called 
VSG consists of weaker and broader features and is likely 
due to PAH clusters (Cesarsky et al. 2000, Rapacioli et 
al. 2005, Rapacioli et al. 2006). The same authors sug- 
gest that PAH clusters could break into free PAHs under 
the effects of a UV radiation field. 

The other th ree observational com ponents were later 
introduced by Uoblin et all ((2008). They were re- 
quired by those authors to fit the spectra of Galac- 
tic and Magallanic Cloud planetary nebulae (PNe) ob- 
served with Spitzer and ISO, which c ould not be fit- 
ted on ly with the templates derived by IRapacioli et all 
(2005). Two of them are broad emission features (BF) at 
8.2/im and 12/zm already seen in some post-AGB stars 
dKwok fc Hrivnaklll989HHrivnak et al.ll2000l: iHony et all 



20011: iPeeters et al.l l2002t Ivan Diedenhoyen et all 120041: 
Szczerba et alll2005l : iKraemer et al.ll2006f) . Peteers et al. 
(2002) classified as C-class objects those in which the 
broad 8.2 /im feature is observed. The 8.2BF and the 
12BF templates were built by Joblin et al. by fitting 
those observed bands in the post-AGB star IRAS 13416- 
6243 MIR spectrum, which is classified as a C-class ob- 
ject prototype by Peeters et al. (2002). It is worth re- 
calling that the contributions of the broad features (i.e. 
8.2BF and 12BF) are important only for the objects that 
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Fig. 6. — Analysis of the PAH spectrum of NGC 6868 after emis- 
sion line subtraction (solid black line) with the template spectra 
displayed in Figure \E\ The best fit is shown as a solid red line, 
while the contribution of each PAH component is shown with the 
same colour code as in Figure \E\ In the legend, we display the 
percentage contribution of the different PAH components to the 
integrated flux in the range between 5.5 and 14 (im. 



are thought to be in transition f r om the AGB to the PN 
phase (e.g. iPeeters et all 120021; Ivan Diedenhoven et all 
12004 iSloan et al.ll2007t lTielens(l2008t Uoblin et al.ll2009l; 
iBerne et al] 120091 ). The second one is a "new" PAH 
type template spectrum, named PAH*, with the 7.8 
/im band shifted to 7.9/im and the 8.55 /im shifted to 
8.65 /im. This te mplate was introduced "ad hoc" by 
iJoblin et al.l (|2008f) . who tentatively associated it with 
large (Nc > 100) PAHs, likely anions, although the large 
cations opt ion could not be excluded (jBauschlicher et all 
l20Ml200flj ). 

With the empirical templates, we obtained best fits like 
that shown in Figure [6] for NGC 6868. Though the best 
fit represents fairly well the general shape of the PAH 
spectrum we find that it underestimates significantly the 
flux of the 11.3 jjum feature. In Figure[7]we show the best 
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Fig. 7. — Fit to the continuum su btracted spectrum of the HII 
galaxy NGC 1482 llSmith et al.H2007l) by using the set of nine dust 
feature template spectra. The meaning of the lines is as in Figure|8] 
Notice that the main contributor to the integrated flux is the PAH 
cations, as is expected in an object with a burst of star formation 
IIAllam andola et al. 1999). The BF contributions are very small. 



fit for NGC 1482. In this case, the MIR spectrum is well 
reproduced with the set of templates, indicating that, 
at odd of our ETGs, the usual PAH spectrum can be 
depicted by a variable combination of these dust com- 
ponents, mainly by PAH+, PAH and VSG. Since the 
7.7/11.3 /im ratio is one of the main peculiarities of the 
observed spectra and the problem affects all our fits, it 
may indicate that the adopted templates, in particular 
the PAH component, are not fully adequate. 

Tielens (2008) points out that the PAH template 
shows more emission in the 6-9 /im region than that ex- 
pected from theoretical calculations and laboratory stud- 
ies of neutral PAHs. The same author suggests two dif- 
ferent possibilities to explain this discrepancy: (a) the 
PAH template could represent the particular mixture 
of PAHs, which could be dominated by neutral PAHs 
but with some, likely small, contribution from ionized 
PAHs, and/or (b) it could also reflect the effect of the 
size on the intrinsic spectrum of neutral PAHs. The lat- 
ter point might indicate that the emitting PAH mixture 
in our ETGs has a size distribution more biased toward 
large and neutral PAHs, than that one characterizing the 
PAH component. 

To bypass these difficulties we also include in our set 
of templates the theoretical spectra of a mixture of ten 
large, compact, and highly symmetric neutral, cationic 
and anionic, PAHs , L-PAH , L-PAH+, L-P AH", respec- 
tively, derived bv lBauschlicher et all (|2008l see their Fig- 
ure 6). In order to compare those theoretical absorption 
PAH spectra with observed PAH emission spectra, band 
shapes, widths and shifts inherent t o the emission pro- 
cesse s must be taken into account (jBauschlicher et all 
120081 120091) . Therefore, w e followed th e recommenda- 
tions of iBauschlicher etaT] ([2008, 2009) concerning the 
fitting of astronomical PAH spectra, and we red-shifted 
the computed spectra by a value of 15 cm -1 , and used a 
bandwidth of 30 cm -1 for the bands short-ward of 9 /im, 
and 10 cm -1 for the bands long-ward of 10 /im. 
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Fig. 8. — Best fits (solid red line) to the PAH spectra of our ETGs (solid black line). The different templates are shown in colour: 
PAH°(solid blue line), 8.2BF component (solid green line), 12BF component (dashed cyan line), VSG (solid magenta line), and Large 
neutral PAHs (solid orange line). In the legend of each panel, we report the percentage of the relative contribution of the different dust 
feature components to the integrated flux in the range 5.5 to 14 fim. Arrows in the spectra of NGC 5044 and NGC 6868 indicate the 
possible detections of the 6.85/im and 7.25/xm features usually ascribed to aliphatic material (e.g. [Chiar "et~aTl 2000: Sloa n et al.ll2007Tl . 



TABLE 5 

Relative contribution in % of the different dust features components to the integrated fluxes of the ETGs. 



Galaxy 


PAH PAH+ 


8.2BF a 12BF a 


BF a 


VSG L 








(8.2+12) 


PAH 


NGC 1297 


61 


9 (11) 26 (30) 


35 (41) 


4 


NGC 5044 


43 5 


10 (11) 26 (34) 


36 (45) 


9 7 


NGC 6868 


35 


12 (13) 36 (38) 


48 (51) 


7 10 


NGC 7079 


58 


8 (10) 24 (29) 


32 (39) 


10 



Notes: 



Values in parenthesis correspond to the relative contribution in % of the broad features components to the integrated fluxes of 
the ETGs in the case of the Virgo stellar continuum template subtraction, see text for details. 



The results of the new fits are shown in Figure [51 where 
the 11.3/im feature is now well-fitted in all cases. In 
Table El we report the corresponding fractional contri- 
butions of the single PAH templates to the best fits. 
The PAH spectra of our ETGs are well reproduced with 
only a few dust feature components: the PAH and L- 
PAH°, the broad feature components (i.e. 8.2BF + 12 
BF) and the VSGs. Particulary stricking is the evident 
deficit of the PAH + emission component, which is the 
main contributor in HII regions and normal galaxies, 



see Figure [7] (e.g. DL07, Tielens 2008 and references 
therein). Even if the PAH and the L-PAH compo- 
nents dominate the total PAH emission, suggesting that 
the main PAH contributors could be likely the neutral 
PAHs, we cannot completely discard the presence of ion- 
ized PAHs in the PAH mixture. That is because, on the 
one hand, the PAH component can contain some con- 
tribution from ionized PAHs (Tielens 2008), and on the 
other, we are using a L-PAH theoretical spectra which 
was computed without including the temperature cas- 
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TABLE 6 

Observed lines ratios, UV magnitudes and derived quantities 





N1297 N5044 N6868 N7079 


[SIII]18.7/[SIII]33.5 
H 2 S(3)/H 2 S(1) 
[NeIII]15.5/[NcII]12.8 
[SIII]33.5/[SiII]34.8 


0.37 ± 0.10 0.35 ±0.15 0.36 ±0.10 0.38 ±0.12 
0.66 ±0.09 1.49 ±0.30 2.85 ± 0.21 1.32 ±0.15 
0.61±0.02 0.52±0.01 0.74±0.01 0.91±0.02 
0.22 ±0.02 0.21 ±0.01 0.80 ± 0.04 0.98 ± 0.07 


m 


20.05±0.12 19.97±0.10 21.10±0.14 


n e (cm 3 ) 

Gpuv(Habing units)'"' 
GFuvT^/no^ 1 / 2 ^ 3 ) 


0.1 - 400 0.1 - 400 0.01 - 400 0.1 - 400 
0.51 0.54 0.20 
< 100 < 110 < 50 



Notes: a The flux of the Habing fields, G FUV = 1, equal to 2.3 X 10~ 3 erg cm" 3 s" 1 at A ~ 1530 A UMathis et al.|[l983T) 



cade effects. iBoersma et af.l (|2010l ) show that the inclu- 
sion of the PAH temperature cascade in the calculations 
affects the relative band intensities. In that case, the 
intensities of the C-H and C-C-C bands increase signi- 
ficatively with respect to the C-C bands (see Figure 18 in 
Boersma et al. 2010), and therefore the fractional contri- 
bution of large ionized PAHs to the PAH mixture in our 
galaxies could be greater than expected from the pre- 
vious analysis. However, even when some contribution 
from ionized PAHs could be expected, the main results 
of our analysis still remains, i.e. the PAH spectra of our 
ETGs are dominated by emissio n from large and neu- 
tral PAH molecules. Moreover, iKaneda et all (|2007al) 
reported unusual 7.7/11.3 /im inter-band ratio in addi- 
tion to no-significant PAH emission at 3.3 /im in the near 
to mid-IR spectrum of the ETG NGC 1316. As the 3.3 
/im emission is usually attributed to emission from small, 
neutral PAHs, Kaneda et al. conclude that the overall 
near to mid-IR PAH spectrum can be explained by emis- 
sion from large and neutral PAH molecules. 

4.2. The Broad Feature components 

The other important contribution that we find in 
our galaxies is the broad feature (BF) components (i.e. 
8.2BF + 12BF) which is, in all cases, higher than 30 
per cent, and even higher than the neutral component in 
NGC 6868. 

Those components are needed to fill the broad 7-9 
/im and 11-12 /um features. In particular, the 8.2BF 
component would fill out the emission excess in the 6 
- 9 /im spectral region mentioned in Section 3.2.1. At 
first glance, it could seem that by subtracting the pas- 
sive ETG template of, we have artificially enhanced the 
broad 7-9 /im feature, and that, as a consequence, we 
need to explain its strength by invoking a broad emission 
feature, identified as the 8.2BF component. First of all, 
our objects are ETGs and a large fraction of their stellar 
populations must be old and quite metal rich (see Table 
IA7l in the Appendix), thus the presence of the old stel- 
lar template is well justified. In addition we adopted a 
template of the passive galaxies in the field, which have 
a less pronounced 8/im dip than Virgo cluster ellipti- 
cal galaxies (see Figure [lj. We also checked how the 
selection of a particular template of an old stellar pop- 
ulation influenced our results. We repeat the analysis 
by using a template based on passive ETG members of 
the Virgo cluster. The results of the analysis are indi- 
cated in parenthesis in Table [5] Finally, the BF compo- 
nent has an even larger broad contribution around 12/im 



(12BF component) and, although in this region the fit 
is not perfect, the presence of this component strongly 
reduces the discrepancy with the observed spectra. As a 
further check we analyzed with the same technique the 
NGC 1482 spectrum (see Figure [7]) which, at odds with 
our ETGs, is well reproduced by a mixture dominated by 
PAH cations, as expected in an HII galaxy (eg. DL07), 
while the BF component is very small. In particular, the 
8.2BF component is negligible. We then suggest that the 
presence of the BF component is not an artifact of the 
analysis. 

So far, the BF features have been detected only in a few 
objects, likely C-stars in transition f rom the AGB phase 
to the PN phase (e.g iBuss et all Il993t IPeeters et al.1 
2001 ISloari et all [20071 : Uoblin et all [2008h . However, 
Sloan et al] (|2005l ) report broad and prominent features 
peaking at ~ 8.3 /xm in the MIR spectra of some 
proto-planetary disks around T-Tauri stars, and sug- 
gest that those broad features are the same as those 
found in the MIR spectrum of e volved stars (C-class ob - 
ject s from Peeters et a l. 2002). iBouwman et al.l (|2008l ) 
and iBerne et al.l (|2009h instead suggest that the broad 
features in T-Tauri stars present different peaks and 
widths than those of the C-class objects. Berne et al. 
(2009) argue that, since the chemical conditions in proto- 
planetary disk and evolved stars are different, there is 
no reason to suppose that the species causing the broad 
features in proto-planetary disks are the same as those 
producing the broad features in evolved stars. Since the 
spectra of our galaxies are well fitted with the parameters 
of 8.2BF and 12BF templates, built from the spectrum 
of the C-class object prototype, we assume that the BF 
emission in our galaxies is caused by the same material 
present in the C-class objects. 

The actual carriers of the BF emission are still un- 
known. These features are commonly attributed to car- 
riers rich in aliphatic material (e.g . iGoto et"ail 120031 : 
ISloan et al.l 120071: Uoblin et alj [20081) . However, C ami 
and coworkers find that the MIR spectrum of the C- 
class prototype object, IRAS 13416-6243 can be well 
reproduced with a mixture of small PAH spectra from 
the NASA Ames PAH spectral database (Cami et al. in 
preparation, see Figure 11 in Tielens 2008). Regardless of 
the actual carrier, which is out of the scope of this work, 
the key point here is that this material seems to be just 
created in the circumstellar envelopes of C-rich stars and 
later processed and/or destroyed by the stronger ISM 
UV fields. This will be later discussed deeper in Section 
6.3. In any case, our tentative detections of the aliphatic 
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bands at 6.85 /nm an d at 7.25 (im (e.g. IChiar et"al"1 
l200?j ISloan <a~ffll200l in the spectra of NGC 5044 and 
NGC 6868 (see Figure [7]) favour the aliphatic material 
interpretation. In that scenario, the BF spectra should 
arise from mixtures of aromatic and aliphatic carbon that 
have not yet been processed by intense ultraviolet radi- 
ation, all owing the more fragile aliphatic bonds to sur- 
vive f e.g.lGoto et al . 2003; S loan et al.ll2007l;lJoblin et al.l 
120081: lTielensll2008|) . A possible evolutionary scenario for 
such features is also proposed: in cool and less evolved C- 
stars (in the post-AGB phase) , the emission is dominated 
by this material. In a hotter ISM environment the fea- 
ture disappears because the aliphatic material has been 
processed into ar o matic mater i al (e. g. iGoto et ail 120031 : 
ISloan et al.l l2007t Uoblin et ai1l2008D . In this way, the 
presence of the BF features in our spectra would indi- 
cate the presence of pristine material. 

5. THE INTERSTELLAR RADIATION FIELD 

The analysis performed above suggests that the un- 
usual PAH inter-band ratios seen in ETGs are due to 
changes in the size and charge distributions of the PAH 
mixture, which would be composed mainly of large and 
neutral PAH molecules. In this section, we consider the 
physical conditions in the ISM and whether or not they 
can account for such a PAH mixture. 

5.1. The PAH charge distribution 
The charge state of PAHs is mainly determined by 

the ionization para meter, GFuvT g as/n c ( see I Bakes et al.l 
120011 : iTiele ns 2005) where Gfuv is the integrated far ul- 
traviolet (900 - 2000 A) radiation field in Habing radi- 
ation field units, T gas is the gas temperature in Kelvins 
and n e is the electron density in cm" 3 . In order to esti- 
mate the values of the ionization parameter for our ob- 
jects, we use the values of the [SIII] 18.7/ [SIII]33.5, the 
H 2 S{2>)/H 2 S{1) line ratios and the GALEX far ultravi- 
olet (A = 1516 A, AA = 268 A) magnitud es within the 
central 5 arcsec from iMarino et al.l (|2010t l , all reported 
in the Table El to quantify the electronic density, the gas 
temperature and the integrated far ultraviolet radiations 
field, respectively. The derived [SIII]18.7/[SIII]33.5 line 
ratios are compatible with mean interstella r electronic 
densi ti es, n fi between 0.1 and 300 cm -3 (Dal e et all 
[200(1 . iTurner et al.l (fl97l showed that H 2 S(3)/ H 2 S(l) 
line ratios ranging from 0.6 to 3 correspond to gas tem- 
peratures of 150 - 400 K . Fi nally, by assuming the dis- 
tances reported in Table IA7[ the FUV magnitudes cor- 
respond to average intensities of the interstellar FUV ra- 
diation field, in Habing units, of Gfuv = 0.51, 0.54, 0.20 
for NGC 5044, NGC 6868, and NGC 7079, respectively, 
which are about half the Galactic value. In the last row of 
Table [6] we also display the corresponding ionization pa- 

1/2 

rameters, GFuvTgks /n e , calculated for those three galax- 
ies. The derived values are in agreement with the ex- 
pected values from the quantification of the 6.2/11.3 fim 
ratio and the 7.7/11.3 /jm r atio as a func t ion of the ion- 
iza tion parameter given by Fla gev et al.l ((20061 ) and/or 
bv lGalliano eTall ((20081) . 

We rou ghly estimate t he fraction of neutral PAHs by 
following (TielensJ ((20051 ) . In this case, the values of 
the averaged ionization parameter of our objects (i.e. 



ture where even the 99% of the PAHs as large as ~ 1000 
C-atoms are neutral. Lower values of the ionization pa- 
rameter would indicate that even larger PAHs would be 
neutral. 

5.2. The PAH size distribution 



GFUvTgas/lle 



100 ) are compatible with a PAH mix- 



iPanuzzo et al.l ((2010) show that a diagram of the 
[NeIII]15.5/[Neir]12.8 vs. the [SIII]33.5/[SiII]34.8 ratio 
allows the different excitation mechanisms in galaxies, 
i.e. starburst photo-ionization, AGN photo-ionization, 
and shock heating, to be distinguished. Using this dia- 
gram, the values of those line ratios for our ETGs, listed 
also in Table [6l are accoun t ed for by the shock+precursor 
models from en et al.l ((2008t ). with shock velocities 

< 150 km s" 1 , pre-shock gas densities between 0.01- 
500 cm -3 and solar abundances, indicating that rela- 
tively fast shocks in a low density mediu m is the primary 
powerin g mechanism in those LINERs. iMicelotta et al.l 
(2010b) analyze PAH processing in shocks by collisions 
with ions and/or electrons, and show that PAHs with 
Nc < 200 are completely destroyed by shocks with ve- 
locities ~ 100 — 150 km s _1 , and only PAHs with 
sizes Nc > 200 could partially survive for a time of 

< 1 — 4 x 10 8 yr. Fast shocks could be originated from jet- 
driven flows or accretion onto a mass ive black-hole, (e.g. 
iNagar et al1l2005HDopita et al.lll997t ). or from turbulent 
motion of gas clouds, either injected by stellar mass loss 
(|Bregman fc Parriottl l2009h or accreted, within the po- 
tential well of the galaxy. However, independently of the 
origin, their presence and the low average UV radiation 
field in the post-shocked interstellar medium could ac- 
count for the changes in the size and charge distribution 
of the PAH mixture in our ETGs. The short lifetimes 
of the PAH molecules in such environments and the fact 
that we see PAH features in the MIR spectra, require 
either mechanisms to protect the PAHs from destruction 
by shocks or a continuous supply of carbonaceous mate- 
rial into the ISM. 

6. DISCUSSION 

Although reminiscent of star forming galaxies, the 
PAH spectra seen in our ETGs show remarkable dif- 
ferences to such objects (see items in Section 3.2). In 
addition, PAH molecules are expected to be destroyed 
in a few 10 6 ~ 7 years in the harsh environment of the 
ETGs. Hence, the presence of the PAHs in the ellip- 
tical galaxies seems to be incompatible with their pas- 
sive interstellar en vironments dominated by hot plasma 
and /o r s hocks (e.g.lBregman et al.l (l992: Mic elotta et al.l 
l2010al |bl: iPanuzzo et al.ll2010t ). So. the key questions to 
answer are (1) why do we see PAH emission in ETGs, i.e. 
what is the source of the PAHs in our objects, and (2) 
why do the PAH spectra present that unusual shape?. 

There are several possibilities to explain this anoma- 
lous emission. We consider them in the following sec- 
tions. 

6.1. Unusual emission due to the presence of an AGN 

ISmith et al.l (|2007l ) noticed that all galaxies of the 
SINGs sample with peculiar 7.7/11.3 fim inter-band ra- 
tios are LINERS, and we also verified that all these galax- 
ies are ETGs. They suggested that the peculiar ratios 
may be due to the presence of low level AGN activity. 
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In that case, the hard UV radiation field would preferen- 
tially destroy small PAH molecules, emitting at shorter 
wavelengths, and leave unaffected large PAHs, which are 
thought to emit mainly at A > 11/im. This effect could 
change the PAH size distribu tion even at large distances 
from the center (jVoitl [l992). Even though this expla- 
nation may be correct to some extent, it cannot be the 
only reason for the observed shape of the continuum sub- 
tracted spectra. A significant emission at 6-9 /im could 
still arise if large PAHs exist in the ionized state (e.g. 
iBauschlicher et al.ll2008l 120091 ; iBoersma et al.|[2010l) . In 
addition, the source of PAH molecules is unknown in the 
Smith et al. scenario. 



6.2. Emitting material recently accreted 

Another p o ssibilit y has been advanced by 
iKaneda et all (|2007bl . I2008t l. They suggest that 
the unusual 7.7/11.3 /im inter-band ratio is due to a 
recent encounter with a gas-rich galaxy that supplied 
the ISM with a dust mixture biased toward large, 
neutral PAHs. A large fraction of neutral PAHs could 
survive the sputtering effects of the hot ISM because 
they were originally formed on the mantle of the dust 
grains. Among their arguments advanced to support 
this possibility we recall: a) the observed missing anti- 
correlation between equivalent width (EW) of PAHs and 
the X-Ray to B-band luminosity ratio, Lx/Lb, implying 
a lack of relation between the PAH sources (stellar mass 
loss rates) and the PAH sinks (X-rays), b) The ob- 
served missing correlation between the 6/mi continuum 
emission and the PAH (11.3/im) emission, implying a 
lack of direct relation between PAH emission and stellar 
populations of the galaxy, c) The correlation between 
PAH emission and continuum emission around 35//m, 
implying a relation between PAH and dust emission. 
Concerning the first two points, we notice that only 
carbon-rich stars (C-stars) can provide carbonaceous 
material, and that these stars are typical of intermediate 
age populations. Their presence in an old system like 
our ETGs would be related to a rejuvenation episode 
and not to the entire stellar population. One would 
therefore not expect a correlation between PAH emission 
and the 6/mi continuum. Furthermore, the EW of PAHs 
do not depend only on the abundance of carriers but 
also on the intensity of the exciting radiation, which 
is mainly UV in the case of the of neutral PAHs (e.g. 
lTielendl200l . This UV flux can v ary from galaxy t o 
galaxy as evidenced by GALEX fe.g. iBuson et al.|[2009D . 
Finally, we notice that point c) only indicates that there 
is a common origin for the exciting flux of PAHs and 
dust. 

Another major difficulty faced by this model is that, in 
many cases, there is no evidence that the last accretion 
episode was sufficiently recent to allow the carriers to sur- 
vive destruction by sputtering or shattering mechanisms. 
Thus, the authors invoke a mechanism that preserves the 
PAHs in the mantle of dust grains. Apart from these con- 
siderations, we notice that it remains unexplained why 
this mixture should be biased toward neutral PAHs, con- 
sidering that the accreted material should be biased to- 
ward the ionized PAHs, as is the ISM in normal star 
forming galaxies. 



6.3. Emitting material continuously supplied by a 
population of C-stars 

The BF component, up until now, has been ob- 
served mainly in carbon-rich pre-plan etary nebulae (e.g. 
iPeeters et atll2002t ISkTan et al.ll2007t) . The prototype is 
IRAS 13416-6243, a carbon star in transition from the 
AGB phase to the planetary nebula phase. This raises 
the possibility that a population of such stars is respon- 
sible of the BF features observed in ETGs. However, the 
comparison with our isochrones shows that this cannot 
be the case because the AGB-PN transition phase is very 
brief, about two orders of magnitude less than that re- 
quired to produce the observed features in the integrated 
spectrum. 

Nevertheless, we advance the hypothesis that a popu- 
lation of C-stars is the source of freshly synthesized car- 
bonaceous material. In our scheme, amorphous carbon 
crea ted within the circumstcllar envelopes of C-rich stars 
(e.g. IWhittetl[l99 2) . is continuously ejected into the ISM 
and then shocked and exposed to a weak UV radiation 
field, as described in Section 5. 

Such a scenario explains in a more natural way the 
origin of the emission in these ETGs. In our case, car- 
bonaceous material is continuously supplied by the mass- 
losing carbon stars, whose origin must be a recent reju- 
venation episode. The timescale for the process to occur 
is provided by the age of the carbon star progenitors, 
which can vary from a few 10 8 yrs to a few 10 9 yrs, the 
latter timescale depending on the chemical composition. 
We thus do not need to advocate a mechanism to pre- 
serve these molecules from destruction (for PAHs several 
10 7 - 10 8 yrs), as Kaneda et al. (2008) did, to explain 
the delay between the actual observation and the epoch 
of most recent merging or of the gas capture episode. In 
our model this delay is simply provided by the stellar 
evolution clock. 

It is interesting to notice that such a population of car- 
bon stars should leave a signature in the MIR spectrum 
of the galaxy, due to t he emission of their carbon rich cir- 
cumstellar envelopes ()Bressan etmil998h This emission 
is similar in shape and peak temperature to the hottest 
residual continuous component left after the subtraction 
of the stellar population template. If this is the case, 
we can make a rough estimate of the number of carbon 
stars responsible for this emission. Assuming for such 
stars L=10 4 L Q and that the light is almost completely 
reprocessed by an envelope of optical depth Ti Mm ~ 10, 
we obtain for NGC 5044 - 15000 C stars. This num- 
ber would require a rejuvenation episode involving ~1 
per cent of the sampled mass of the galaxy, about 1 Gyr 
ago, which is a typica l value found among field galaxies 
(jLonghetti et alJl2000f ). 

In summary, our model can explain both the low 
7.7/11.3 /im ratio and the presence of BF families as well 
as relax the strong constraint imposed by the sputtering 
time on the epoch of the gas accretion. 

Thus we argue that the anomalous PAH spectrum is 
a direct signature of the presence of an intermediate age 
carbon star population. 

7. SUMMARY AND CONCLUSIONS 

We have analyzed the MIR spectra of four ETGs char- 
acterized by an unusual PAH emission spectrum. These 
galaxies, which belong to the sample that will be fully 
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presented in lPanuzzo et al.l (|2010D . show prominent PAH 
complexes at 11.3 /im and at 17 /im, with abnormally low 
7.7/11.3 /im inter-band ratios. 

In contrast to star forming galaxies, the IRS spec- 
tra of our ETGs are dominated by stellar emission at 
short wavelengths. Thus, in order to investigate their 
continuum subtracted spectra, we have first subtracted 
realistic templates of old stellar populatio ns, derived by 
high s ignal-to-noise Spitzer observations (Bressan ct al. 
I2006al) . After passive stellar template subtraction, we 
have performed a detailed analysis of the individual dust 
emission components. 

A surprising result is that the 7-9 /im spectral re- 
gion requires features which are different from those 
found in normal late- type galaxies, such as those dis- 
cussed by DL07. In this spectral region some important 
features look significantly wider and/or shifted toward 
longer wavelengths. 

With the aim of understanding all the above peculiar- 
ities we have analyzed each spectrum by mea ns of dust 
feature templates (see e.g. iJoblin et al.l |2008|) and the - 
oretical PAH spectra (see e.g. iBauschlicher' etani200l . 
This method has proved to be very efficient for the under- 
standing of the nature and origin of the PAH emission in 
Galactic objects ,but, to our knowledge, has never been 
used for external galaxies. 

The result of this analysis shows that the ETG spec- 
tra are dominated by two main components: (1) the large 
and neutral PAH component, responsible for the peculiar 
7.7/11.3 /im and 6.2/17 /im ratios, and (2) the so called 
BF component, which produces the broadening and the 
shift of the features in the 6-9 /im spectral range. The 
PAH cations, that are the dominant family in normal 
galaxies, are almost completely absent in these spectra, 
indicating that the PAH emitting mixtures in our galax- 
ies have anomalous size and charge distributions, biased 
to large and neutral PAHs. In our sample, the BF com- 
ponents contribute approximately 30-50 per cent of the 
total PAH flux between 6 and 14 /im and are responsible 
for the observed shift and widening of the fitted Drudc 
profiles with respect to DL07. To test whether this re- 
sult is an artifact of the adopted procedure, we have also 
analy zed the SINGs HII galaxy NGC 1482 (jSmith et al.1 
2007). In this case, the MIR spectrum is well repro- 
duced by a mixture dominated by cationized PAHs, as is 
expected in an HII galaxy, while the 8.2BF contribution 
is negligible. 

Up to now, the BF features have been observed mainly 
in evolved carbon stars and have be en associated with 
pristine carbonaceous material (e.g. Honv et all 
Peeters et all 120021 : iBerne et al.l 120071 



Sloan et all 



2001 



2007 



Joblin et al.ll2008l ). Thus, our analysis provides convinc- 



ing evidence that we are seeing pristine carbonaceous 
components, still relatively unprocessed by the ISM en- 
vironment. The most natural explanation we can ad- 
vance is that this material originates from a population 
of carbon stars. 

We stress that in the above analysis, abandoning the 
usual black body template for the underlying old stellar 
populations, in favour of a more realistic one derived from 



passive ETG galaxies, was a fundamental step. 

In order to account for the anomalous PAH mixture, 
we characterized the interstellar radiation fields in our 
objects by using mid-IR line ratios sensitive to the elec- 
tron density, gas temperature and excitation mechanism. 
The results of this analysis suggest that shock models, 
with velocities between 100 - 150 km s~ 1 , in a low den- 
sity ISM, can account for the line ratios observed in our 
galaxies. If such shocks are present in our objects, they 
would first destroy the smaller PAHs on timescales of a 
few 10 8 yrs, and would help to maintain the observed size 
distribution. The extremely low values of the ionization 
parameter estimated for our objects lead to a predomi- 
nantly neutral charge distribution. 

Present data do not allow us to identify the actual 
impact of shocks on the dust size distribution, but the 
peculiar properties of the PAH emission in ETGs may be 
the result of an environment where the PAH molecules 
are processed in shocks, and excited by the weak UV 
radiation field of the old stellar population. 

We propose that the unusual PAH spectrum arises 
from the combination of two effects. First, a rejuvenation 
episode that recently ceased gives rise to a population of 
carbon stars that is now continuously feeding the ISM 
with fresh carbonaceous material. From the comparison 
of the residual warm dust component seen in the galaxies 
with the circumstellar envelope dust emission of a typi- 
cal carbon star, we estimate that the rejuvenation frac- 
tion is about 1 per cent of the sampled mass. We notice 
that kinematical studies of NGC 5044, NGC 6868, and 
NGC 7079 show star/gas counter-rotation; a strong indi- 
cation of a past accretion event (see appendix). Second, 
shocks process the ejecta of these carbon stars, destroy- 
ing smaller PAHs. The weak UV radiation field of an old 
stellar population excites the processed material, main- 
tains a suitable ratio between neutral and ionized PAHs 
and results in the anomalous spectral features. 

It remains to be clarified whether the presence of these 
shocks and the occurrence of the rejuvenation episode 
are causally linked. While there is a general consensus 
that AGN and star f ormatio n activity are linked (see e.g. 
iBressan et al.l 120021 (2006b), the situation is much less 
clear for LINERS. However, if our interpretation is cor- 
rect, our results indicate that LINERS could also have 
hosted star formation activity in the recent past, suggest- 
ing a LINER/post-starburst connection. 
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APPENDIX 



RELEVANT PROPERTIES OF THE GALAXIES SO FAR 

The four galaxies considered in this paper belong to the sample of 40 ETGs in low density environments investigated 
in the mid-infrared by Panuzzo et al. ( 2010). This sample is biased towards ETGs having emiss ion lines in their optical 
spectra (see the analysis produced bv lRampazzo et aL1l2005t [Annibali et al.ll2006l 120071 [2010K It is known, however 
that t he presence of ion ized gas in ETGs is quite common both in the field as well as in clusters (see e.g. lMacchetto et al.l 
fl99& IS arzi et~aT][2005l) . Most of these ETGs are classified as LINERs using optical diagnostic diagrams (see Annibali 
et al. 2010). 

In the Lick-IDS line-strength indices analysis, performed by I Annibali et all (|2007f ). the four ETGs discussed in this 
paper were assigned quite different ages. 

NGC 1297 is the oldest with a luminosity weighted age of 15.5±1.2 Gyr, but only scant information exists in the 
li terature about t his ga laxy. 

lAnnibali et al.l (120071) attributed to NGC 5044 the very uncertain age of 14.2±10 Gyr. The galaxy, located in a 



rich group of galaxies (|Tullvill988[ ). is rich in dust in the central 10" with a clu mpy distribution. The ionized gas has 



a filamentary structure with an extension of about 40" (|Macchetto et al.lll996l ). The gas velocity profile is irregular, 
with many humps and dips, while the inn er (within 1/3 of the effective radius) stellar velocity profile is counter- 
rotating with respect to the outer r egions (ICaon et all |2000|) . This galaxy, then, is a very peculiar object, being a 



possible merger/accretion remnant. iRickes et al.l ( 20041 ) studied the ionized gas component in NGC 5044 suggesting 



the presence of both a non-thermal ionization source in the central region and an additional ionization source (possibly 
hot post-AGB stars) in the outer parts. 

NGC 5044 has a set of infrared observatio ns. It has been detec ted by IRAS and by ISO (|Ferrari et alj|2()()2. and 
reference therein) . The galaxy belongs to the iKaneda et al.l (120081) sample and t hey detected PAH features and also 
H2 rotational emission lines and ionized species using Spitzer~ ]Temi et all (|2007f ) discussed IRAC and MIPS Spitzer 
observations of NGC 5044. They report interstellar emission at 8/j.m detected out to 5 kpc as well as 70/im emission 
from cold dust exceeding that expected from stellar mass loss. In view of the short sputtering lifetime for extended 
dust (w 10 7 yr) they conclude that the extended dust cannot result from a recent merger with a gas-rich galaxy. They 
support the view that the complex and highly fragmented dust clouds are highly transient and created by stellar mass- 
loss in the central «1 kpc. The dust clouds are intermittently disrupted and heated by energy released by accretion 
o nto the central black h ole (AGN Feedback). 

lAnnibali et~aT1 (|2007l) provide luminosity weighted ages for NGC 6868 and NGC 7079 being 9.2±1.8 Gyr and 6.7±1.1 
Gyr, respectively. 

The Fabry-Perot observations of NGC 6868 ((Pi ana et all 19981) show that t he line-of-sight velocity field of the ionized 
gas component has a velocity amplitude of ± 150 km s _1 . ICaon et alJ (|2000l ) show that along the axes at P. A. =30° and 
70° the gas and stars h ave similar kinematica l properties, but along P.A.=120° the gas counter-rotates with respect to 
the stellar component. Zeilinge r et al.l ([1996) noticed the presence of an additional inner gas component which they 
suggested could be due to the superposition of two unresolved counterrotating components, one dominating the inner 
region, the other dominating the outer parts. Stars also show a kinematically-decoupled counterrotating core. The 
stellar velocity dispersion decreases towards the galaxy center. The above kinematics suggest that NGC 6868 could 
h ave had a recent accretion episode which could explain the relatively low luminosity weighted age. 

iBettoni fc Gallettal (fl997l ) found that in NGC 7079 the gas is rotating in a direction opposite to that of the stars 
(gas counter-rotation) . 

We may conclude that, according to the current interpretation of the peculiar gas vs. stellar kinematics found in 
NGC 5044, NGC 6868 and NGC 7079, accretion/merging e vents could have occurred i n the se galaxies. 

Both NGC 6868 and NGC 7079 has been observed in CO. iHuchtmeier fc Tammannl (|1992f ) observed with the SEST 
NGC 6868 and provide a CO upper limit. Assuming this value we obtain an upper li mit to the cold H2 m ass of 
2.0 x 10 8 Mq (adopting a gala ctic X factor X co = 4.6 M Q (K km s _1 pc 2 ) -1 , see ISolomon et al.l fl987h . For 
NGC 7079. IBettoni et alJ (|2001h estimated a mass of cold H 2 of 8.82 x 10 7 M within the central 2 kpc radius. 

In Table [ATI we collect relevant properties of the galaxies. Column (1) gives the name of the galaxy; column (2) 
the morphological type; column (3) the class of activity; columns (4), (5) and (6) the lumin osity weighted age and 
the metallicity and a-cnhancement computed from optical spectra bv lAnnibali et afl ([2007) . Column (7) lists the 
molecular masses obtai ned from CO observati ons and by using the galactic X-factor. Column (8) displays the X-ray 
luminosity reported in (O' Sullivan et al.ll2001l) . In column (9) we give the luminosity distance taken from NED using 
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TABLE A7 

General properties of the sample 



NGC Type Activity Age Z [at/Fe] M C H2 log L x D 

Class (Gyr) (M ) (erg s- 1 ) (Mpc) 

1297 SABO pec: LINER 15.5±1.2 0.012±0.001 0.29±0.04 19.8 
5044 E0 LINER 14.2±10.0 0.015±0.022 0.34±0.17 42.74 42.8 

6868 E2 LINER 9.2±1.8 0.033±0.006 0.19±0.03 < 2.0 x 10 s 41.23 37.5 

7079 SBO LINER 6.7±1.1 0.016±0.003 0.21±0.05 8.82 x 10 7 34.2 

Notes: (Col. 2) Morphological types arc derived from RC 3; (col. 3) Activity Class from Annibali ct al. (2010); col.s 4-6 from 

lAnnibali et~aU (2007); (col. 7) H 2 content from Huch tmeier k. Tammannl fl99^ 1 and lBettoni et al.1 1 120011) for NGC 6868 and NGC 7079 
respectively; (col. 8) X-ray Luminosity from IOgIe^tani|20*0*6T h (Col. 9) Distances adopted are derived from NED. 



H =73 km s^Mpc" 1 . 



